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In order to study further the promising free piston Stirling engine architecture, there is a need of an ana¬ 
lytical thermodynamic model which could be used in a dynamical analysis for preliminary design. To aim 
at more realistic values, the models have to take into account the heat losses and irreversibilities on the 
engine. An analytical model which encompasses the critical flaws of the regenerator and furthermore the 
heat exchangers effectivenesses has been developed. This model has been validated using the whole 
range of the experimental data available from the General Motor GPU-3 Stirling engine prototype. The 
effects of the technological and operating parameters on Stirling engine performance have been investi¬ 
gated. In addition to the regenerator influence, the effect of the cooler effectiveness is underlined. 

© 2010 Elsevier Ltd. All rights reserved. 


1. Introduction 

The Stirling engine, invented in 1816, running according to a 
reversible closed cycle knew a practical use as a reliable and sure en¬ 
gine, during almost one century before being supplanted by the 
spark-ignition engine. Nowadays, the Stirling machines are in com¬ 
mercial use only as heat pump, used mainly for cryogenic cooling 
and air liquefaction. As an engine, the Stirling remains a field of 
numerous researches and development works. Recent experimental 
realizations demonstrate power densities as well as noteworthy effi¬ 
ciencies [1]. One of the most promising applications of the Stirling 
cycle is the free piston Stirling engine (FPSE) configuration [2-4]. 

The optimal design of FPSEs is a difficult task. Indeed, no 
mechanical linkage fixes the strokes and phase angle for the mov¬ 
ing elements. Hence, a global dynamic analysis is required to pre¬ 
dict the periodic steady operation. Due to the complexity of this 
analysis, the isothermal assumption is usually adopted. Hence, 
the pressure as a function of the piston and displacer positions 
can be expressed in an analytical way. Linearization methods are 
then used to obtain the performances of the engine [5-8]. How¬ 
ever, these models do not take into account the thermal losses of 
the engine which lead to erroneous predictions of the perfor¬ 
mances. Therefore, there is a need for an accurate analytical ther¬ 
modynamic isothermal model which can be used in accordance 
with the dynamical analysis of the FPSE for preliminary design 
purpose. 
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Many investigators have studied the effect of some heat losses 
and irreversibilities on the engine performance indices. In the 
many parameters to be taken into account, dead volume as well 
as non-ideal regeneration have the highest influence on the Stirling 
performances compare to all the technological parameters of a 
practical engine [9,10]. Popescu et al. [11] showed that the most 
significant reduction in performance is due to the non-adiabatic 
regenerator. Kongtragool [12] studied the influence of the regener¬ 
ator efficiency and the dead volumes on the work as well as the 
efficiency of the machine. However, this study does not include 
the heat transfers through the temperature difference at the heat 
source and sink. 

On the basis of the conventional entropy techniques, for the 
studying of solar Stirling engine cycle performance, Costea et al. 
[13] included the effects of heat transfers, incomplete heat regen¬ 
eration and irreversibilities of the cycle as conduction, pressure 
losses or mechanical friction between the moving parts. Timoumi 
et al. [14] developed a precise second order model which includes 
all the losses at the same time. The method based on a lumped 
analysis approach leads to a numerical model and has been used 
for the optimization of the General Motors GPU-3 [15,16]. 

But this type of models appears to be not suitable for a prelimin¬ 
ary design stage of FPSE. Nevertheless, the extensive study based on 
the GPU- is a reference model to validate new developments. 

Stirling machine with dead volume can be analytically studied 
using the Schmidt approach [17] as far as isothermal evolutions 
and ideal regeneration are assumed. Following a second order ap¬ 
proach according to Martini [15] classification, the Schmidt results 
can be completed by an energetic balance including a non-ideal 
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Fig. 1. General schemes of a Stirling machine with thermal interfaces. 


regeneration. In addition to the previous limitations, the thermal 
transfers between the machine and its surrounding thermal 
sources have to be taken into account in a Stirling machine model. 
They induce great thermal losses which have to be limited in tech¬ 
nological realizations. Heat exchangers at the hot and cold ends 
transfer the thermal energy to the working fluid. The necessary 
thermal gradient to induce this transfer gives way to a temperature 
difference between the thermal source and the working fluid at the 
considered part of the machine. The Newton’s law allows the mod¬ 
elling of the relation between the thermal transfer rate and the 
temperature difference. In the work of Senft [18], a constant effec¬ 
tiveness of thermal exchangers with respect to fluid flow rate is as¬ 
sumed. Given that on the one hand, the heat transfers are related to 
the operating frequency and on the other hand the frequency ap¬ 
pears to be a main optimization parameter, the study of Senft must 
be examined in more details. Consequently, the model developed 
in this study accounts for the effect of the flow rate to evaluate 
the thermal effectiveness of the heat exchangers. 

The present work consists in the elaboration and the validation 
of an analytical isothermal model as a first step to the development 
of an accurate global analysis of FPSE. The obtained analytical mod¬ 
el can be used to analyse Stirling engines including their dead vol¬ 
ume, imperfect regenerator and external in addition to internal 
thermal transfers. 

Firstly, the analytical Schmidt results are recalled. With the 
addition of the thermal flaws, energetic characteristics of a Stirling 
engine are set out. Then, the heat exchangers and the regenerator 
efficiencies as a function of fluid flow rates are added. Though it is 
not a free piston type engine, the GPU-3, it is one of the most doc¬ 
umented Stirling engine and its data are therefore used to compare 
the model results to the experimental ones. For a wide range of dif¬ 
ferent pressures and frequencies, the model shows a good correla¬ 
tion with the experimental data. The model is also compared to the 
classical adiabatic second order results and shows good agree¬ 
ments with them. In the last part, an optimization is performed. 

2. Thermodynamic analysis 

In this study, the conventional schematic representation of a 
“thermodynamic machine” (framed part in Fig. la) is completed 
by three thermal coefficients. Denoted a, b and c, they are associ¬ 
ated with the heat exchangers efficiencies for the first two as well 
as conduction phenomenon for the latter. 

The Newton’s law for expressing heat losses will be adopted to 
model the thermal transfers in the machine. 

2A. Stirling engine characteristics 

The generic schemes of a Stirling machine as well as the tem¬ 
perature distribution are described in Fig. 1. The displacer takes 
place in the left end whereas piston is at the opposite side. 


In order to obtain an analytical expression of the engine pres¬ 
sure, the usual following assumptions are adopted: 

• Working fluid within the engine chambers stays at constant and 
unvarying temperatures. 

• The temperature of the working fluid within heat exchangers 
volumes V h and V k are Tu and T L respectively. 

• Temperature within the regenerator can be described by a linear 
evolution between T a and T L . 

• The regenerator behaviour is symmetrical. 

• The ideal gas law is adopted. 

• The pressure is the same throughout the machine for each con¬ 
sidered position of both the piston and the displacer. 

• Harmonic movements for the piston and the displacer. 

Organ has demonstrated that an equivalent Stirling machine al¬ 
ways exists [19]. As a result, an equivalent Stirling machine what¬ 
ever alpha, beta or gamma Stirling engine mechanical arrangement 
considered can be assed using the same one-dimensional geometry 
model. Therefore, we choose here to develop our model on such an 
equivalent Stirling machine for which expansion and compression 
volumes evolutions are: 

V e = Vd/2{\ + cos 0(f)) 

V c = kV d /2{\ + cos(</>(t) - a)) 1 j 

In which 0(t) is the angle position and a the phase angle. The piston 
and displacer swept volume ratio is denoted k = V P IV d . 

The total mass of gas inside the various chambers is defined as 
the sum of the mass of the gas contained in each of the chamber 
defined in Fig. 2: 

m = m HC -hm h -hm R -h m,- + m k + m cc (2) 

If we substitute the expressions of the masses in (2) with the ideal 
gas law (e.g. m HC = we obtain: 
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Fig. 2. Permissible area for L 
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In which R is the considered gas constant per unit of mass and T R the 
mean temperature of the regenerator. 

The Schmidt analysis results [17] are recalled hereafter: 


mR 1 

s 1 +/?cos(0(t) - 0) 



with: 


• Thermal efficiency 

The losses associated to the regenerator have the highest influ¬ 
ence on the performances. We anticipate that the shuttle losses as 
well as the gas spring hysteresis have little influence on the perfor¬ 
mances. By doing this, these losses are neglected compared to the 
former ones. 

The thermal efficiency of the engine is defined by the ratio of 
the available power by the added thermal power: 


P = 


V^T/CCOSa + T 2 + K 2 


tan 9 


2v + K + T 

/csina 
/ccosa + T 

1 


(5) 


1 = 2 ^t_ 

s V d 2V + /C + T 

_ Vcc VrTl Vhc 

“ Vd" Vd" 

In which t = TjJTu. 

Consequently, the main characteristic parameters of the engine 
are given by the following analytic expressions: 


V 


Mean effective pressure 
mR 1 


mean 


V 1 -/* 2 

Heat added (Qe) and heat rejected (Qc) for each cycle: 


( 6 ) 


Qe= pdV HC =p mean 

J cycle 


TCV d 
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sin0 


(7) 
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[ pdV C c = -Pmean-JT- 

J cycle P 
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1 - 
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1-J8 3 


sin(a - 0) 


( 8 ) 




Qe stands for the heat added to the “hot part” of the engine and Qc 
for the “cold part”. 

• Efficiency and power 

The ideal thermodynamic efficiency is exactly the Carnot 
efficiency: 


^ = 1 - T l /Tu 


(9) 


Useful mechanical work can be determined from Eqs. (7) and (8) 
and varies notably as a function of the mean effective pressure: 


Iyl/ 17 uk(t- l)sina (, [\ Tf 

Wi — Pmean^d ,—~— =>——~- ( 1 \/l ft 

PVK 2 + T 2 + 2/crcosa 


( 10 ) 


Thus, the indicated power can be expressed as the product of the 
mechanical work by the operating frequency /: 


Pi =jw, 


( 11 ) 


2.2. Thermal approach 


Pth 


Pth 

Q.h + Or + Q.k 



In which Qt is the conduction loss, and Q R is the thermal power re¬ 
lated to the regenerator inefficiency. 

We define the various thermal impedances ratios: 


Pcond ~ C cond/ a ( 14 ) 

p R = m R C v /a (15) 


where rh R is the fluid mass rate inside the regenerator. 
Hence: 


Qr = ap cond Tu(t - T) (16) 

Qr = (1 - ejapJuP - t) (17) 


The regenerator effectiveness e depends on the temperatures of the 

Jf T 

engine. We define e = T'u < Tu the highest actual outlet 

temperature of the fluid from the regenerator. It must be noted that 
as far as a symmetrical regenerator behaviour is assumed a single 
effectiveness is defined and e can be conversely defined as 
e = with T[> T L the actual lowest temperature of the fluid 
from the regenerator. 

Finally the thermal efficiency is expressed as a function of the 
non-dimensional parameters: 


Pth 


1 + Sy - { — 

1 -£-(Pcom/ + (l -e)P R )« 1 -T) 



• Thermodynamic conditions 

The second law of thermodynamics requires that the thermal 
efficiency does not exceed the Carnot efficiency that is: 

J7 t „<l-T (19) 

For a machine which operates as an engine the available power P th 
is always positive. Then we set: 

Pth > 0 ( 20 ) 

Therefore, using the two previous conditions (19) and (20), the per¬ 
missible values of £ are within an hatched domain shown on Fig. 2. 

In order to obtain the maximum performances an optimal case 
is chosen. Thus, the inequality (19) can be switched to equality. 
Consequently, the optimal ratio £ op tim between the heat source 
temperature and the temperature of the expansion chamber of 
the engine can be obtained. Using Eq. (18): 


• Thermal power 

The cycle average power turns out to be given in an alternative 
way. According to the power balance of the machine as described 
as in Fig. la, the cycle average power can expressed as: 

P* = aT H (l+8r-i-Zfo) (12) 

In which £ = T V /T H is the temperature ratio between the heat source 
and the highest temperature of the engine, 3 = b/a and r = T C /T H . 


* _ sr -\-t 

Qoptim — 2 v Z1 7 

(<5 + 1)t - (p cond + (1 -e)p R )( 1 - t ) 

2.2.2. Expression of the thermal coefficients 

According to the Newton’s law the coefficients a and b are as¬ 
sumed to be representative of the heat exchange for the hot and 
the cold end respectively. Thermal effectiveness of the heat 
exchangers are dependant on the machine parameters and more 
specifically the operating frequency /. Therefore, it is important 
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to consider this dependency to deepen to the work of Senft [18] 
where a constant thermal coefficients are assumed. 

The a and b thermal coefficients are deduced from a simple con¬ 
vection model relations: 

q = h k A wk 

b = h k A wk { } 


In which h h and h k are the convection heat transfers, A wh and A wk are 
the wetted areas of the heater. 

Based on experimental correlations, h is expressed as a function 
of the Colburn J-factor J h [20]. 


h(Re) 


hPr 213 

cm 



In which Pr is the Prandtl number and Af f the free flow area and the 
cooler respectively. C p is the constant pressure specific heat capacity 
of the working fluid. 

The evolution of J h with respect to Re is of the form 
Jh(Re) = Ci Re~ n . For tubular heat exchangers which are used for 
the GPU-3 Stirling engine, the Colburn J-factor for heat transfer is: 


Re ^ 3000 J h = exp(0.337 - 0.812 log(fie)) 

3000 ^ Re ^4000 J h = 0.0021 

4000 ^ Re ^ 7000 J h = exp(13.31 + 0.861 log(Re)) 

7000 ^ Re ^ 10000 J h = 0.0034 

10000 ^ Re J h = exp(-3.575 - 0.229 log(Re)) 


A polynomial interpolation for J h will be used in the model. 

Let us now move to the Reynolds number Re. The mass flow rate 
can be estimated from de Schmidt analysis and the Organ works 
[19]. Fig. 3 plots the Reynolds number for the heater and the cooler 
as a function of x for the GPU-3 parameters. The usually reported 
turbulent behaviour of flow is predicted. 

The thermal coefficients can eventually be evaluated from the 
engine geometric and operating parameters. 


?m — 


eo — 1 — 


NTU 
2+NTU 
1 


?DM 


NTU 

NTU 


as proposed by Martini [15 

is suggested by Organ [19] 

2 FR 


f(FR ) where/(FR) = 1 - exp(-NTU) 


2 +NTU NTU (2 +NTU)' 
if FR ^ 1 and f(FR) = 1 - exp(-NTU/FP) otherwise, 
by De Monte [21] 


(24) 


The NTU can be evaluated through the correlation of the Colburn 
factor for a mesh grid regenerator with the Reynolds number [15]. 
Once again, the calculation of the Reynolds number for the regener¬ 
ator relies on a coupled Schmidt - Organ approach. As shown in 
Fig. 4, laminar flows occur in the regenerator for the considered 
operating frequencies. 

The results for the various expressions of e with respect to the 
temperature ratio x and for three operating frequencies are plotted 
in Fig. 5. 

The Organ definition gives more important values than the two 
others. The flush ratio has little influence on the regenerator effi¬ 
ciency. As a consequence, the simple Martini expression will be 
used in the following. 


• Conduction losses 


Most of the thermal conduction appears through the regenera¬ 
tor which separates the hot and cold areas. We choose an expres¬ 
sion of the heat flux related to the conduction losses based again on 
the Newton’s law and a constant coefficient c C0nd defined as: 

^ _ b M Affreg 1 - % ni:\ 

L cond — I'wi'conn j m 

^reg 1 y 

In which N conn is the apparent connectivity related to the thermal 
connexion in the mesh grid regenerator [19] and k w the thermal 
conduction of the regenerator grid material. 


2.2.2. Expression of losses 


Mechanical effectiveness 


• Reheat loss 

A non-ideal regenerator can not provide an outlet fluid temper¬ 
ature equals to the adjacent chamber temperature. It follows that 
the hot heat exchanger must provide an extra heat called Qr de¬ 
fined by Eq. (17). 

Simplified expressions of e with respect to NTU can be found in 
the literature. They all involve NTU and some include additional 
operating parameter such as the flush ratio FR [21]. 

The Eqs. (24) hereafter give different expressions of e: 


Save for thermal losses, mechanical losses must also be taken 
into account to evaluate the brake output power and global 
efficiency. 

On the one hand friction occurs within kinematic linkages joints 
for the piston and displacer as well as for the output shaft, on the 
other hand moving parts seals can lead to great friction forces. Both 
count in a single mechanical effectiveness coefficient called p mec . 
Moreover we assume a fixed value of r\ mec whatever the operating 
conditions. 


/= 20 Hz 



Fig. 3. Evolution of the Reynolds numbers for the heat exchangers with respect to 
temperature ratio t for three operating frequencies. 
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Fig. 4. Evolution of the Reynolds number for the regenerator with respect to 
temperature ratio t for three operating frequencies. 
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Fig. 5. Regenerator efficiency e with respect to temperature ratio t for three 
operating frequencies. 

Another source of mechanical dissipation lies in pressure drops 
associated to the fluid flow trough exchangers and regenerator de¬ 
noted A p exch and A p R respectively. Pressure drops can be calculated 
using the friction factor coefficient related to the Reynolds number. 
Therefore, the total mechanical effectiveness is: 

Pdiss = Imec + A Pexch/P + ^Pr/P (26) 

Senft has demonstrated the drastic effect of the buffer pressure on 
the mechanical effectiveness [22]. He established the central theo¬ 
rem as will be recalled hereafter. 

The additional mechanical features are to be considered as shown 
on the schematic representation of Fig. 6. W out is the brake output 
power and arrows stand for the power transfers between parts. 

One can finally give the expression of the indicated power: 


P j — W i — W exp W comp 

(27) 

And the output power: 


Pout = w 0ut = T] diss (W + - W-) 

(28) 


The global final cyclic mechanical effectiveness f] MEC is the ratio of 
the output power P out to the indicated power P,. 

The central theorem established by Senft is recalled hereafter: 


The cyclic mechanical effectiveness of any engine with volume 
compression ratio r = V max / V min , temperature ratio z = T L / T a 
and mechanical effectiveness s< E, cannot exceed v\ M ec- 

1 -E 2 

tjMEc(FT,r) = E -S(T,r) 

S(T,r) = 0i/ir^ 1 

_ T lOg(T) - ( 1 + T) (log( 1 + t) - log( 1 + r) - log(r)) , 

( ’ (1 - T)log(r) 7 

(29) 


This upper bound is the mechanical efficiency of an ideal Stirling 

engine buffered at its optimum constant pressure and having a 

mechanism of constant effectiveness E. 

Substituting rj diss for E in the above theorem, the cyclic mechan¬ 
ical efficiency rj MEC can be added to the previous model. Fig. 7 plots 
the mechanical efficiency with respect to t. For the high tempera¬ 
ture ratio of considered here, the effect of the buffer space pressure 
can be neglected. If low temperature differences are considered 
which is the case solar application the mechanical efficiency curve 
shows a sharp slope downward. 

3. Results and discussions 

The operating point can be obtained from the equality of indi¬ 
cated and thermal power given in Eqs. (11) and (12) respectively. 
Fig. 8 plots parametric curves of non-dimensional powers 
Pi = Pi/(fPmean V d) and P* h = P th /(fp m ean V d) for various operating fre¬ 
quencies. The curves of the two powers with respect to tempera¬ 
ture ratio z intersect at one point. 

Because of the enhanced heat exchange effectiveness, P* th in¬ 
creases as/increases and the operating point moves to lower val¬ 
ues of z. However, if / is greater than 90 Hz, the thermal transfer 
properties reach a limit resulting in an almost constant apex for 
P* h . The conclusions are quite different from the constant heat 
transfers case studied by Senft, in which the temperature ratio z in¬ 
creases with / in order to balance heat transfer governed by the 
Newton’s law. Therefore, the temperature differences between 
external and internal temperatures always increase. 

3.1. Model validation 

The model results are compared with those obtained by Urieli 
and Berchowitz [9] and Timoumi [14] for the same conditions of 
the GPU-3 (i.e. working fluid Helium, mean pressure 
Pmean = 4.13 MPa, and operating frequency/= 41.72 Hz). The re¬ 
sults presented in Table 1 shows a good agreement with adiabatic 
models which take into account the pressure drop, the regenerator 
efficiency and conduction losses. 

As expected from the isothermal assumption, the predicted effi¬ 
ciency is 13% higher and the power and the work by cycle are very 
close when compared to the adiabatic model without shuttle and 
gas spring hysteresis losses. 

Most of the studies deal with a single operating point for the val¬ 
idation stage. Because the model is dedicated at design and optimi¬ 
zation, it is compared hereafter to GPU-3 experimental values for a 
wide range of parameters especially for the pressure and frequency. 

The geometrical parameters as well as operating data for Gen¬ 
eral Motor GPU-3 are summarized in Table 2. 



Fig. 6. General Stirling machine scheme with thermal interfaces buffer pressure and mechanical parts. 
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Vmec\% ] /= 45 Hz 

- /= 60 Hz 



/= 5 Hz 



Fig. 8. Evolution of the operating point for various frequencies. 


To allow easier comparison, the performances of the model are 
expressed in original units from the GPU-3 experimental test re¬ 
sults [15]. The Fig. 9 shows the variation of the torque, brake out¬ 
put power and efficiency with respect to operating frequency for 
various mean pressures. The dotted lines are from the experimen¬ 
tal results whereas the straight lines stand for the model results. 


Table 1 

Results from different models. 


Type of model 

Heat 
input, Jl 
cycle 

Indicated 

output 

W 

power 

Jl cycle 

Thermal 

efficiency 

(%) 

Adiabatic model 

327 

8286.7 

198.62 

62.06 

Urieli and Berchowitz [9] 
adiabatic model 

— 

8300 

— 

62.5 

Timoumi dynamic model 
without losses [14] 

314 

7109.3 

170.4 

54.96 

Urieli and Berchowitz quasi¬ 
steady flow 

— 

7400 

— 

53.1 

Timoumi dynamic model with 
loss dissipation (Ml) 

291 

6372.4 

152.47 

53.3 

Urieli and Berchowitz quasi¬ 
steady flow (pressure drop 
included) 


6700 


52.5 

Timoumi dynamic model Ml 
and internal conduction loss 
(M2) 

294 

6355.2 

152.32 

52.64 

Developed analytical isothermal 
model 

258 

6087 

146 

52.9 

Timoumi (M2) and external 
conduction loss 

314 

6061 

145.27 

46.94 

Timoumi dynamic best model 

262 

4273 

99.5 

38.49 

Experiment 

- 

3958 

- 

35 


The model shows good agreement with the whole data. The 
analytical model qualitatively represents the trends of the curves. 
Due to the isothermal assumption, the efficiency from the model 
is above the experimental one as expected from this type of models 
[15,23]. The efficiency results show a constant discrepancy of 
about 46% with the experimental ones. 

The experimental results show a static torque which prevents 
the engine from running for the lowest value of the mean pressure. 
This effect is not taken into account in the model. Consequently, 
the torque results are forced to zero by subtracting the value ob¬ 
tained for the lowest mean pressure. By doing this for the torque 
and the power, a constant difference of about 25% for each fre- 
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Fig. 9. Brake output power and efficiency curves with respect to operating 
frequency and pressure. 


Table 2 

Parameter values of the GPU-3 Stirling engine. 


Heat source T H (°C) 

900 

Cold source T c (°C) 

15 

Phase angle a (°) 

120 

Swept volume ratio k 
(-) 

Working fluid 

1.01 

Swept volume V d (cm 3 ) 

120.88 

Hydrogen 

Mean pressure p mean 

0.69 < p < 6.89 

Mechanical 

90 

(MPa) 

Operating frequency / 

20 </< 50 

effectiveness p mec (%) 


(Hz) 

Wire mesh 


Regenerator length 

22.6 

regenerator 


(mm) 


Wire diameter d w (pm) 

40 

Porosity of the 
regenerator matrix % 

0.759 

Regenerator hydraulic 

0.03 

Material conductivity 

16.6 

radius r hR (mm) 


(Witt 1 K 1 ) 
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(b) P[kW] 



Fig. 10. Influence of the fluid mass brake output power and efficiency. 


P [kW] 



Fig. 11. Influence of the mean pressure on the brake output power - efficiency 
curve. 

quency is obtained. It is worthy of note that the reduction of the 
torques values by increasing the operating frequency is well repro¬ 
duced by the model. 

3.2. Optimization 

3.2A. Effect of the fluid mass 

The total mass of the gas is related to the mean pressure by the 
Eq. (6). The proportional relation between the power and the mean 
pressure can be seen in Figs. 10a and 11. 


This result is close to the results obtained for the GPU-3 and the 
effect of the fluid mass on the efficiency is different for the work of 
Timoumi [14]. An optimal fluid mass of m = 1.6 g (i.e. a mean pres¬ 
sure Pmean = 6.89 MPa) for which a maximal value of 54% efficiency 
can be inferred (see Fig. 12b) which is 2% higher than for the pro¬ 
totype reference value shown in dashed line in Fig. 10. In the same 
time, the power is raised from 6.09 kW to 10 kW. 

3.2.2. Combined effect of fluid mass and frequency 

For a given mass of gas, the operating frequency can be opti¬ 
mized to improve the performances. Fig. 13 shows the evolution 
of the power-efficiency curves for various pressures. The GPU-3 
reference operating value of 41.72 Hz appears to be close to the 
optimal value whatever the pressure. 

3.2.3. Effect of the regenerator length 

By varying the regenerator length, the performances are modi¬ 
fied accordingly. The additional dead volume as well as an increase 
of the pressure loss make the power reaches a maximal value. 
However, for small regenerator length, i.e. below 15 mm, the 
power as well as efficiency quickly decreases, as shown in Fig. 12 a. 

The power-efficiency curve shown in Fig. 12b underlines the 
choice of the optimal power optimization chosen for the GPU-3 
operating parameters. 

3.2.4. Regenerator thermal conductivity 

The conductivity of the material constituting the regenerator 
matrix has strong effects. Fig. 13 shows that with an increase of 
the matrix regenerator thermal conductivity leads to a reduction 
of the performance due to an increase of internal conduction losses 
in the regenerator. For a value close to the steel thermal conductiv¬ 
ity (/< w ^40 Wm" 1 K _1 ), the optimal efficiency is about 19% at a 


(a) p [kW] n[%] (b) p[kW] 




Fig. 12. Influence of the regenerator length on the brake output power and efficiency. 
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Fig. 13. Influence of the conduction loss on the brake output power - efficiency curve. 


P [kW] 


140 Hz 



Fig. 14. Influence of the cold heat exchanger performance on the brake output power - efficiency curve. 


higher frequency than for the stainless steel used for the GPU-3. 
For a low conductivity material such as the Inconel 625® 
(k w « 9.8 Wm" 1 I< -1 ), the efficiency can be increased to 57%. 

3.2.5. Effect of the cooler efficiency 

Although water cooling is used for the GPU-3 engine, a Stirling 
engine can use free or forced air flow for the cooler. Therefore for 
technological design, one can be interested in cold heat transfer 
coefficient appraisal (see Fig. 14). 

An increase of 10% of S leads to a sensitive variation of both the 
efficiency for a given frequency. Conversely, a decrease of 10% in¬ 
duces the same absolute difference on the efficiency. Again, the 
major role of heat transfers in Stirling engine is outlined here. 

4. Conclusion 

As a first step for design and optimization of FPSE, an analytical 
model has been elaborated. To deal with the discrepancy between 
the high theoretical efficiency of Stirling and a constructed proto¬ 
types the major losses are be taken into account. Moreover, they 
are related to the geometrical and physical parameters of the pro¬ 
totype design. 

In order to obtain a simplified analytical model, an isothermal 
assumption is adopted. The model integrates the regenerator effi¬ 
ciency and conduction losses, the pressure drops and the addi¬ 
tional heat exchangers effectiveness. Thus, the influence of the 
operating parameters on the global engine performance can be as¬ 
sessed. As a validation stage, we applied the parameters of the 
GPU-3 engine data on the developed model. Despite using the iso¬ 


thermal assumption and simplified heat transfer model, the results 
are close to the experimental data and are in agreement with clas¬ 
sical models. 

An optimization of these parameters has been carried out using 
the GPU-3 engine data. It is shown that a reduction of losses and a 
notable improvement in the engine performance can be reached. In 
addition to the well known influence of the regenerator, the drastic 
effect of the cooler effectiveness is underlined. 

As a conclusion, a handy model is developed and can be used in 
a first design procedure and optimization with respect to a specific 
application. 
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